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Abstract 

We estabhsh the relation between two objects: an integrable system related to 
Painleve II equation, and the symplectic invariants of a certain plane curve T^tw de- 
scribing the average eigenvalue density of a random hermitian matrix spectrum near a 
hard edge (a bound for its maximal eigenvalue). This explains directly how the Tracy- 
Widow law Fgue; governing the distribution of the maximal eigenvalue in hermitian 
random matrices, can also be recovered from symplectic invariants. 
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Motivation 



The famous Tracy- Widom law [32], describes the statistics of the largest eigenvalue of 
a random hermitian matrix of the GUE ensemble [28] . 
Consider the large deviation function 

Fn{,o) = lnProb[Amax < a] for a < 2 (0-1) 

Here 2 = limAr^oo E(Amax)- It is proved, that J^n{,o) has an expansion in l/N"^, either 
in [31] by probabilistic methods, or in [16] by Riemann-Hilbert asymptotic analysis for 
orthogonal polynomials [5, 13, 14, 15]: 

J-^ (a) = + J] Ar2-2s (a) (0-2) 

with F^{a) independent of N. It is known that the F^{a) are the symplectic invari- 
ants introduced in [18] (see Section 1.1 below). To be short, they can be computed 
recursively, purely by means of algebraic geometry on a certain plane curve S(a). 

In [7], we have computed those coefficients F^{a), also for the genealization to /3 
ensembles (/3 > 0) with polynomial potential. In this article, we shall consider only 
the hermitian case (/3 = 2). Because of universality, we may restrict ourselves to a 
quadratic potential, i.e. to the GUE ensemble. 

Then, heuristically, one can study not so large deviations by pluging a = 2 + N~'^/^s 
{s < 0) in Eqn. 0-2. The result is that (see [7]), with this scaling, each term of the 
expansion becomes of order 1: 

FN{a = 2 + ^ Aou^ + ^(-5/2)3(1-5) (Stw) (0-3) 

Ague is a constant that one can compute, and F^(Stw) are the symplectic invariants 
of a certain plane curve Stw? which is a limit of the curve S(a) when a — )■ 2. This 
limit object is in fact universal (it is the same if we start from a generic potential), it is 
the curve of equation = x + 1 — 2, up to normalization of x and y. Heuristically, it is 
expected that the RHS of Eqn. 0-3 is the exact s — )■ — oo asymptotic of limAr_>.oo J^Ar(a = 
2 + iV~^/^s). Indeed, this method for /3 > reproduces all previously known results 
[11, 30, 12, 2] on the left tail asymptotics of {3 Tracy- Widom laws. It also predicts the 
constant for any /3 > 0, and provides an algorithm to compute the expansion at all 
orders. In our case, for GUE: 

PauE(^) = 2V-eC'(-i) exp (-^ - ^ + Y.{-s I2f^'-^^ F^S^w)) (0-4) 

\ S>2 / 
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From Tracy and Widom [32], we know that PqueI-^) = Fgue(s), where Fgue is related 
to Painleve II equation, which itself appears in relation to an integrable system. How- 
ever, by performing a expansion and computing F^, the integrability structure 
is not manifest. At first, it was not clear why the asymptotic series in the RHS of 
Eqn. 0-4 should be related to Painleve II. In this article, we fill this gap and give a 
proof of the following: 

Proposition 0.1 

-^ + E(-^/2)'^"^^^^(Stw)) (0-5) 

is the asymptotic when s — )■ — oo of a Tau function for two 2x2 compatible systems of 
ODE: 

s) = L(a;, s), ds^{x,s) = M.{x,s)'^{x,s) (0-6) 

whose compatibility condition is equivalent to a Painleve II equation: 

q"{s) = 2g3(s) + sq{s) (0-7) 

The identity of Eqn. 0-3 was based on a heuristic argument, the double scaling limit 
of a matrix model, which would require a mathematical proof. In this article, we start 
from the RHS given by heuristics (i.e. from the world of the topological recursion), 
and we prove that it is indeed the expansion of the LHS. We do not use a Fredholm 
determinant representation of Prob[Amax < o]- 
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Outline 



In the first part of tlie article, we gather older facts relevant for our problem: 

• the framework of the topological recursion and its properties (Section 1.1-1.2). 

• the general relation between 2x2 integrable systems and loop equations (Sec- 
tion 1.3-1.5). 

In a second part, we refine the technical hypothesis underlying this very relation 
(the outcome is Theorem 2.1, our main theoretical result). We explain what really has 
to be checked (Section 2.1), and we sketch how this be done in practice (Section 2.3). 

The third part is devoted to the application of these ideas to an integrable system 
related to Painleve II (Sections 3.1-3.2). Subsequently, we prove Prop. 0.1 (Sections 3.3- 
3.5). 

Let us mention that many ideas used in the proof also appeared in the work of 
M. Cafasso and O. Marchal [9] with different purposes (they related the (2m, 1) minimal 
models appearing in the merging of two cuts in the 1 hermitian matrix model, to a 
hierarchy of equations containing Painleve II), and in fact have their origin in the works 
of the second author with M. Bergere [3, 4]. In this article, our goal is to close (in the 
hermitian case only) the alternative approach to Tracy- Widom laws presented in [7]. 
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1 Basics around the topological recursion 



1.1 Definitions 

The topological recursion and symplectic invariants were axiomatically defined in [18], 
and we refer to [19] for a review, and to Appendix B for a summarized definition. 
It consists in an algorithm associating some numbers and differential forms 

uj^(T,){zi, . . . , Zn) to a regular spectral curve S. 

For our purposes, we call spectral curve the data of: 

• a plane curve (C, x, y), in other words a Riemann surface C, with two meromorphic 
functions x and y, C C 

• a maximal open domain f/ C C on which x is a coordinate patch, called physical 
sheet. 

• a Bergman kernel B{zi, Z2), i.e. a differential form in zi & C and in Z2 G C, such 
that, in any local coordinate ^: 

B(z.,z.) = #f^ + 0(l) (1-1) 

The zeroes of dx are called branchpoints (name them Oj), and a spectral curve is 
said to be regular when these zeroes are simple and are not zeroes of dy. In other 
words, when a/x — x{ai) is a good coordinate around G C, and when y behaves like 
y{z) ~ y{ai) + y'{ai) ^J x{z) - x{ai) with y'{ai) ^ 0. 

We shall not give the full definitions of F^(E) = Co'o(^) '^n(^) here, and rather 
refer to Appendix B or [18]. Their construction is axiomatic and relies only on algebraic 
geometry on C. The construction is made by recursion on 2g — 2 + n. 

We rather mention the essential properties that we use here: 

• F^(S) is invariant under any transformation (x, y) — )■ (xo, yo) such that |dxAd?/| = 
|dXo A dyo\. For this reason, the F^(Il) are called symplectic invariants. 

• For 2 - 2^ - ri < 0, ^^(S) G T*(C) ® . . . ® T*(C), i.e. wS(S)(zi, . . . , z„) is a 
meromorphic differential form in each Zi G C, symmetric in all ZiS, and it has 
poles only at the branchpoints, of maximal degree 2(3(7 + ^ — 2), with vanishing 
residue. 

• for 2 — 2g — n < 0, w^^(S) are homogeneous of degree 2 — 2g — n, i.e. if we change 
y \y, we have a;^(S) ^ A^-^f-" u;9(S), and F9(S) ^ A2-29F9(s). 
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• i^^(S) have nice scaling properties when the spectral curve approaches a singular 
spectral curve. We will be more precise when needed. 

• <^^(S) have nice properties under variation of the spectral curve. If is a 1- 
parameter family of spectral curves, let us write {dtydx — dtxdy) — and 
represent the differential form Q with the Bergman kernel: 

n{z)^ f MOB{^:^) (1-2) 

JQ* 

(this is form-cycle duality: Q* is the cycle dual to the differential form ^2, and in 
case the cycle is not regular, it is accompanied by a Jacobian Aq). Then: 

dtUjf,{Et){z,,...,Zn)= [ Af2(0^^+l(St)(e,^l,...,^n). (1-3) 

1.2 Matrix models and topological recursion 

The topological recursion allows to find all order asymptotics of the large N expansion 
of matrix models. Let us recall how it works for the 1 hermitian matrix model with one 
hard edge, relevant for the application to Tracy- Widom law. For any potential V such 
that I J^d.^e"^^^^^] < oo for some positive e, let us consider N real random variables 
Ai, . . . , Ajv with joint probability law of density: 

1 ^ 
d/^iv(A) = — — TT |A,-A,f •TTdA,exp(-ArT/(A,)) (1-4) 

^ ' l<i<j<N i=l 

The probability of having max Aj < a is given by V{a) — ZjsfQ — oo, a])/Zjv(oo), where: 



Zjv(J)= / d/.iv(A) (1-5) 

dfj^N is the probability distribution induced on the eigenvalues of a random hermitian 
matrix of size N xN, with probability law dM exp (— yTr V{M)^ up to normahzation. 
dM here is the canonical Lebesgue measure on the real vector spaces of hermitian 
matrices. The observables invariant under conjugation are spanned by the cumulants: 

N / N \ 
i=n \ji=l ' ■'V 

When lnZ7v(a) admits an asymptotic expansion in powers of we write: 

ZN{a) = C^rexpf J]Ar2-25ir9 J (1.7) 

\9>0 / 

Wn{xi,...,Xn) = J]Ar2-29-W^^^(a;i,...,x„) (1-8) 
9>0 
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where Cat is a normalization constant depending only on A^, and and W-^ are 
independent of but depend on a and V. Before coming to the coefficients themselves, 
let us say that when the density of eigenvalues has a large N limit, p(^)d.^, which admits 
a single interval Jq C] — oo,a[ as support, it is proved that such an expansion exists 
[31, 16]. If the support consists of several intervals, no such expansion exists (there are 
fast oscillatory terms depending on A^), but there is still a way to extract the partition 
function from the topological recursion [20]. Yet, this discussion is out of the scope of 
this article. The theorem proved in [17, 10] is that when the large expansion exists, 
then: 

F9 = F^{E{a)) (1-9) 
W^'\x{zi),...,x{zn))dx{z,)---dx{zn) = c^f (S (a) ) (zi , . . . , z„) (1-10) 

for the spectral curve S(a) determined by the equation y = V'{x)/2 — Wi{x) and the 
Bergman kernel determined by: 

B{zi, Z2) = (w^ix{zi),xiz2)) + , , , ..A dx(zi)dx(2;2) (l-H) 

y{x) is analytic on the x complex plane, except on the eigenvalue support Jq, where it 
admits a discontinuity y{x + iO~^) — y{x — iO~) = 2i7rp(x). C is a Riemann surface on 
which y can be analytically continued. 

1.3 2x2 differential system and topological recursion 

To any 2x2 first order differential system: 

d^^ = L^, t ^) (1-12) 



is associated the integrable kernel: 



Xi - X2 

We shall restrict ourselves to TrL = 0. So, dxidet"^) = 0, and we can choose the 
normalization: 

det ^ = 1 (1-14) 

In [3] were also introduced correlators Wn{xi, . . . ,Xn) and connected correlators 
Wn{xi, . . . ,Xn). The connected correlators are defined by: 



Wi(x) 



lim (lC(x,x') ^ 1 (1-15) 

x'^x V X - X' J 
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yV2{,xi,X2) = -/C(a;i,a;2)/C(a;2,xi) - 7 r^j (1-16) 

{xi - 

n 

>v„(xi,...,xo = (-1)"+' n^(^*'^-«) (1-1'^) 

(T cycles of Sn i=\ 

and the correlators by: 

Wn(xi,...,x„) ="det"/C(xi,x,) (1-18) 

where "det" means that each occurence of )C{xi,Xi) in the determinant should be 
replaced by >Vi(xj), and each occurence of IC{xi,Xj)IC{xj,Xi) by —W2{xi,Xj). In other 
words the >V„ are the cumulants of the Wn- 
For example: 

Wi = ^ 9J - ^9,0 = -{d,^ 4> - dj^ 0) (1-19) 

Eqn. 1-14 implies that all correlators are symmetric in the Xj's. It can be checked that 
they do not have poles at coinciding points i ^ j- The spectral curve of a first 

order differential system is defined by the plane curve C of equation: 

det{y - L(x)) = 0. (1-20) 



Theorem 1.1 (proved in [3]) Assume that: 

(i) L depends on some parameter N, and has a limit when N ^ oo. 

(ii) The spectral curve of the system Eqn. 1-20 has a large N limit which is 
regular, and has genus 0. 

{iii) >V„(xi, . . . , Xn) admits an asymptotic expansion when N ^ oo of the form Wn = 
^^>g A^i-2g yyg^ ^^^^ 2g — 2 + n > 0, >V^(xi, . . . , Xn) may have singularities 
only at branchpoints o/Sqo- 

Then, the expansion coefficients of the correlators are computed by the topological re- 
cursion applied to the spectral curve defined by Sqo with Bergman kernel B{zi,Z2) = 
{W^{x{zi),x{z2)) + l/{xizi) - x{z2)f)dx{zi)dx{z2): 

yV^{x{zi), . . . ,x{Zn))dx{zi) ■ --dxiZn) = U^(J2oo){Zi, . . . , Zn) (1-21) 

We stress that it is the large N limit of the spectral curve (also called "classical", or 
"semiclassical" spectral curve) which is relevant in this theorem. A similar result holds 
when Soo is not of genus under an extra hypothesis, but this will not be needed here. 
So far, we did not speak of an integrable system, the theorem itself tells something 
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about one differential system alone. However, in many examples and in this article, 
proving that (iii) holds is made possible by the existence of other compatible systems 
dt^"^ = Mj\E', where tj is a parameter of L, and is rational in x. So, we shall see 
in Sections 2.1-2 A that hypothesis {Hi) can be considerably weakened. 

1.4 Generalities on the Schrodinger equation 

Let us write: 

L = f ^ M (1-22) 



c —a 

Eqn. 1-12 imply that ip and are two independent solutions of a second order equation: 

dl^p + Bd^%l) + V%1) = {) (1-23) 

where: 

B = (1-24) 

V = detL -d^a + a^ (1-25) 

All the same, ip and are independent solutions of the second order equation: 

dl^ + B d^i) + U i) = Q (1-26) 

where: 

B = (1-27) 

c 

V = detL + d^a-a— (1-28) 

c 

We shall use a classical change of function in the theory of Schrodinger equations to 
study the asymptotics of \E' at singularities of L. We are grateful to M. Bergere who 
shared his experience on that subject with us. We can always write: 

^ = ^ ""^^ {- f J^) ^^'^^^ 

Then, Eqns. 1-23-1-26 are equivalent to: 

\dlf = 2U{dJ) + {d,U)f (1-31) 

^dlJ = 2U{dJ) + {dJJ)f (1-32) 
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with U, U given by —V, —V up to a schwartzian derivative: 



4 2 c 



1.5 Tau function and topological recursion 

Here, we modify slightly the systems by introducing a parameter in front of each 
derivative (very often, turns out to be redundant with other parameters). For a 
large class of 2 x 2 first order compatible systems (for example this includes the (p, 2) 
minimal models which are finite reductions of KdV) of the type: 

\-d^m = l.m \-dt^! = m.m (1-35) 

and under the hypothesis of Thm. 1.1, it was stated in [4, 9] without details that 
exp ^^^>Q A^^~^^ F^(Soo) j is a formal Tau function. We shall make the argument 
explicit here. We are considering any system like Eqn. 1-35, such that L and are 
rational fractions in x. We assume that L has a large N limit L*^°\ and an expansion 
in powers of 1/A^. 

1.5.1 Definition of the r-function 

Here, we consider that is a fixed parameter. To any solution \&(x,t) of Eqn. 1-35, 
Jimbo-Miwa-Ueno [26] associate a Tau function t{1^ as follows. Let us consider the 
behavior of \E'(x, t) near each pole Xa of L(a;, t), and write the x — )■ Xq asymptotics 
(valid in some angular sector near Xq): 

^(x, t) = ^„(x, t) exp (T«(x, t)) (1-36) 

where e^°' contains the essential singularity and the pole of \E'(x,t), and \E'a(x,t) is 
analytical when x — )■ x^. Then, the compatibility of the system Eqn. 1-35 implies that 
there exists a function r(t) such that: 

dtATiT = -VResdxTr^-i(9,^J(9t^Tj (1-37) 

a 

= -V Res dxTr(^-^a^^)e"'^"(9..T^)e^" (1-38) 

f ^ X — ^x^ 



x—^x 

a 



Proposition 1.1 Assume thatTa{x,i) is family of diagonal matrices, and write = 
Ta^s + cl (where = diag(l, —1) ). We have: 

dt^ In r = 2 V Res dx (dtT^) Wi (1-39) 



X — ^X(y 

a 



with the correlator Wi(x) = {dx4')4> — {dxip)(p introduced in Section 1.3. 
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proof: 

^^^9^.^ = ^^^L^ is traceless like L, and Ta{x,i) is a family of commuting ma- 
trices. Thus, only the component of T of the Pauli matrix is relevant. Eqn. 1-38 
yields: 

dt^ Inr = -2 V Res dx (9tT„) (1-40) 

a 

And we compute using det \E' = 1 : 

^ I ^ M d^ij j 



^(^x-^) - ^xV^ -V^ 5x0 + <9x0 
j->V]_ ^ d^(f) - (f) d^(f 



:i-42) 
:i-43) 



□ 



1.5.2 Large limit 



We use Section 1.4 and keep the same notations, except that now, each derivative 
comes with a power of 1/A^. This modification can be easily traced back: 



^P(x) = cteW^^^^expf-AT (1-44) 



where / and / satisfy: 



^ dlf = 2U{dJ) + {d.,U)f (1-46) 



2Ar2 



^ dlf = 2UidJ) + {d.U)f (1-47) 



2Ar2 



with the potential: 



U = _detL+-(c)..-a^j + -(-4^--^) (1-48) 

U = -detL+- -9.a + a^ +— (1-49) 



A^ V c J N"^ \A 2 c ^ 

Let us write with superscript ° the large A^ limit of these quantities. In particular: 

L'-(:;:r;) 

11 



Recall that we defined the spectral curve through the equation det{y — L(x)) = 0. 
Since L is traceless, it is equivalent to y'^ = — det L(x). Hence, the large limit of the 
spectral curve has equation: 

Soo(t) : = -detL(°)(x) (1-51) 

Let us define y{x) by this equation, with a choice of the branch of the square root 
depending on the asymptotic x — )■ we require to define ip (see Eqn. 1-54 below). 
According to Eqns. 1-48, the potential in the large N limit is given by: 

[/(o)(a;) =Z7^°\x) = i/2(a;) (1-52) 

And, from Eqn. 1-46, we get: 

/(°)(a;) cx 1/yix), 7^'\x) cx 1/yix) (1-53) 
Matching the asymptotics gives the proportionality constants and we get: 



^W(x) = cte J^^exp ( -N / yiOd^ ] , (1-54) 



= ^ y exp (^N y(Ode) . (1-55) 
1.5.3 1/A^ expansion of the r-function 

We now want to study asymptotics when N ^ oo. We assume that 

T^{x) = NT^;^\x) = -Na^ (^j^ y^i^ (1-56) 

where (■ ■ ■)-,q, is the divergent part when x — )■ Xq, and o is an arbitrary origin of 
integration. In other words, we assume that L(x) and L*^°^(x) have the same poles 
{xq,}, that the behavior of \E' at these poles is entirely given by the large A^ limit of the 
system. This happens quite often in practice. Then, we rewrite Eqn. 1-39: 

dt^r = -2N ^ Res d^^ym^^ m{x) (1-57) 

Provided the hypothesis of Thm. 1.1 hold, the topological recursion can be applied, 
with the spectral curve (Soo(t)) defined in Eqn. 1-51. Let us compare to J-'(t) defined 
by the large A^ asymptotic series: 

^(t) ^ 5^Ar2-29^.(s^(t)) (1.58) 

9>0 
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We also define: 

W,{x) ^ 5^iVi-2^iyf(S^(t))(x) (1-59) 
<?>o 

Let us call C, the underlying curve of Soo(t), and B its Bergman kernel. We shall use 
the properties of the topological recursion (Section 1.1) to compute the variations of 
J^(t) wrt tj. We have to represent 

{dt^ydxm=ResB{z,0 I ?/dx (1-60) 

Jo 

Notice that the RHS do not depend on o. For 2x2 systems, Soo(t) is hyperelliptic. So, 
each pole x^ of ^ ydx has two preimages and Za in C. One can move the contour 
to surround these points: 

{dt^ydx){0 = -J2 Res_ 5(^,0 / dt^ydx (1-61) 

a ^" 

Subsequently: 

9t^F5(Eoo(t)) = -T" Res_ dx{z)Wf{z) [ dt^ydx (1-62) 

a " 

In the case of an hyperelliptic curve, the involution is a symmetry in the construction 
of Wf, and one finds that the residues at Za and equal. Hence: 

= -S^ Resdx{z)2NWi{x{z)) f dt ydx 
^ Jo ' 

= dtAnr (1-63) 

This proves the identification of the resummation of symplectic invariants with the r 
function. 



2 Refinement on the assumptions 

2.1 A remark on hypothesis {in) 

Starting from any solution to j^dx"^ = L\E', with L a given 2x2 matrix, rational 
in X and traceless, the connected correlators W„ were defined in [3] and Section 1.3 
such that they satisfy "loop equations". When these loop equations are written in 
terms of the r function, they appear identical to Virasoro constraints'^. Loop equations 

■^Independently of [3], Virasoro constraints have also been found in [27] in the case of 2 x 2 
Schlesinger system. 
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can have many solutions, but their solution is unique if one requires some analyticity 
properties. In fact, provided there is a large parameter in the problem, and provided 
an expansion: 



exists, the topological recursion gives the unique solution of the hierarchy of loop 
equations satisfied by W^, whose analytical properties are specified by the data of the 
spectral curve. 

The matter is then, for the one is interested in, to find the appropriate spectral 
curve Soo- It is often easy to prove that has an expansion in and to discuss 

its analytical properties (see Section 2.3). However, hypothesis (iii) of Thm. 1.1 is 
stronger. Indeed, the two essential properties of Eqn. 2-1 are: 

• Fixed parity: only powers of 1/A^ with parity (—1)" appears in the expansion of 



• Factorization property: the leading order of Wn is 0{N'^ "). 

We believe that in practice, "fixing parity" is only a matter of rescaling correctly with 
N the parameters of the differential system^ dx'^ = L\E' and that L*-*^^. In Appendix A, 
we give an example based on a Lax pair associated to Painleve IIq equation, where the 
fixed parity property does not hold if one does not rescale a as well. 

Now, we argue that the factorization property is automatic, as soon as one realizes 
dx"^ = L\E' as the member of a full integrable hierarchy. We present this argument for 
2x2 systems with 1 pole at x = oo (the case of interest in this article), but most of 
these ideas could be adapted without difficulty to several poles and larger systems. We 
are currently working out the generalization of all arguments presented in this article 
to the case of d x d differential systems [8]. 

2.2 Loop insertion operator and properties 

In some sector near z = oo, we consider a solution \E' whose asymptotic is given by 
\E' = ^^0,6^^°"^, with \i/ regular at x = oo, and: 



for some constants tj. It is a standard result in integrable systems [1] that tj can be 
used to define commuting fiows. The result yields a matrix L and solution \E', which 

^This is not a void statement, because Theorem 1.1 asks by convention for a in front of each dx 
in the differential system, and also that lim^v-yoo L exists. 




(2-1) 




(2-2) 
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now depends on all t,-, with sectorwise asymptotic given by Eqn. 2-2, and such that: 



— (9^^ = 
N 

N ^ 



Mj{x) 



(2-3) 



[■■■]_ ^ denotes the divergent part when a; — )■ oo, and obviously is a polynomial in 
X of degree j. Let us define the loop insertion operator as a Laurent series of operators 
when X — oo: 

6. = T^ (2-4) 



^-^ x^~ 

It is merely a way to collect all the fiows^. 

Lemma 2.1 We write ~ the functions with their essential singularity re- 

moved. The following 2x2 matrix: 

1 - $(x) 0-3$"^' 



P x) 



-ip{x)4>{x) ip{x)(t){x) 
-i/j{x)(j){x) ip{x)(t){x) 



2(P(X2)-P(xi)) 

X2 - Xi 

[P(xi),P(a;2)] 

X2 - Xi 



$(xi) 



(2-5) 

(2-6) 
(2-7) 



is a projector, and we have: 

1 

N 
1 

iV 
proof: 

Eqn. 2-5 is a straightforward computation. We just notice, according to the represen- 
tations of Eqns. 1-29-1-30, that the essential singularities cancels when we compute the 
cross-products: ipcf) = ip cf) and ipcl) = ip (p. Then, we observe that P is a rank 1 matrix: 



with: 



-1 

1 



'(0 0)(V^ V^)J 



J^ = -l, 



J-1 



(2- 



(2-9) 



By definition TrP = 1. Hence, P^ = (TrP)P = P. Now, let us compute the action of 
the loop insertion operator on dt^"^ = Mj\E' translates into: 

1 



a.^(xi) = '^{xi)x{a3'^ \xi) ^(a;i) - ^(xi)a;V3 



(2-10) 



5 J_ 

N 



Sx given by Eqn. 2-4 gives a realization of the formal loop insertion operator introduced in [3] 
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It is implied that [• ■ •]_ is the divergent part when xi — )■ oo. Doing the summation 
over j: 

1 



X2 — Xi 



^(Xi) - ^(Xi) 



X2 — X\ 



X2 — Xi 

2(P(X2)-P(xi)) 



^(Xi) 



X2 - Xi 

We apply this result to compute: 



^5,,P(xi) = -^5,, (^(xi)a3^(xi) ^ 



P(X2) -P(X2)7, , , 7, , ,P(x2)-P(xi) 

^-^(Xi)a3^ l(xi) + ^(xi)(T3^ '(Xi) 



X2 - Xi X2 - Xi 

-(P(X2) - P(Xi))P(Xi) + P(X1)(P(X2) - P(Xi)) 
X2 — Xi 

[P(X1),P(X2)] 



X2 - Xi 



□ 



Corollary 2.1 The integrable kernel is self-replicating: 

^ 6x^)C{xi, xs) = -/C(xi, X2)/C(x2, xs) (2-11) 
Proposition 2.1 T/ie correlators can be expressed in terms of the projector P: 

Wi{x) = -ArTrP(x)L(x) (2-12) 



lTr(P(xi) -P(x2)) 



2 



anc? /or n > 3.- 

W4x....,x.) = iV-(-l)-^ ^ ^ TrP(^i)P(x.a))--P(x..-H,) ^ 

(2-14) 

Or concisely with the loop insertion operator: 

Wnixu ...,Xn) = iV-" 6,,--- In r (2-15) 

As a consequence of the existence of a loop insertion operator and Eqn. 2-15, the 
factorization property is automatic, and the requirements of hypothesis {iii) have to 
be checked only on Wi(x). 
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2.3 Analyticity properties of the expansion 
2.3.1 Summary of the problem 

Let us start with two 2x2 matrices L(x,t) and lSA{x,t) such that: 

• L and M are traceless, and rational in x. 

• They are solution of the compatibility equation: 

l^iL - ^dM + [L, M] = (2-16) 

• They admit an asymptotic expansion when — )■ oo of the form: 

L = ^A^-'l(^), M = ^ iV-' M(') (2-17) 

Z>0 l>0 

and consider the equations: 

— = L^, — (9t^ = (2-18) 

N N 

In other words, we first restrict ourselves to a hierarchy with only one time t. We shall 
discuss the analytic properties of ip (the discussion would be similar for ip, and 0). 
From the differential system wrt x, we know (Section 1.5.2) that it admits a large N 
asymptotic expansion of the form: 

^{x) = ( 5^iV-'^«(a;) I e-^/^^^)^? (2-19) 
\i>o J 

where o is some integration constant which may depend on t. The leading order is 
given by: 

,7". ^ ^ (2-20) 

and y is determined by the large spectral curve associated to the differential system 
in x: 

Soo : = -detL(°) (2-21) 

Similarly, we may define the large spectral curve associated to the differential system 
in t: 

; y_' = -det M(°) (2-22) 
We shall see that Eqo and S^^ play an important role to locate the singularities of ip. 
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2.3.2 Analytical properties 

Let us describe the analytical properties of ■?/'*^°-' in the variable x: 

• It has branchpoints (order ^-^^p^, r G Z) at branchpoints of Soo- 

• It has branchpoints (order r G Z) at the zeroes or poles of which are 
not branchpoints of Sqo- 

We can insert Eqn. 3-41 in one or the other differential system to determine recursively 
the subleading orders. It is more convenient to write the two second order equations 
satisfied by ip: 

-^dl^lj+^Bd^^ + Vij = 
^d^ij+^mi^ + Vi^ = 



where: 



b ^ 

V = detL-is.L„ + lL„^ = V + 5:^'-''^"' 

'^^ 1>1 

V = detM-i-9,Mn + ^Mn^ = V + $^iV-'Z« (2-23) 
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1>1 



Knowing we can compute recursively, if we assume that B^^\B^^^ are not identi- 
cally zero: 



= ^[-a,2V'(n(a.i/)^« + 2y(a.^«) 



m=0 



(2-24) 



m=0 



(2-25) 



A priori, ■?/;('+^) has singularities at points where ■?/;(™') (for < m < /) had already 
singularities, and at poles of L and M. Besides, if we consider the differential system 
wrt X, new singularities may appear at each step: 
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• At the branchpoints of (Soo) (through 1/y). 

• At the zeroes (in the x variable) of B'^'^\ 

If we rather consider the differential system wrt t, new singularities may appear only: 

• At the branchpoints of (through 1/y). 

• At the zeroes (in the x variable) of B^'^\ 
Subsequently, il)^^\x) may have singularities only: 

• At the poles in the expansion of L and M. 

• At the branchpoints of Sqo- 

• At the branchpoints of S^^. 

• At the common zeroes (in the x variable) of B*^^^ and B^'^\ 

In many examples (and in this article), the fourth point is the key ingredient to rule 
out singularities at the zeroes of y which are not branchpoints. It is often enough to 
have one time t, but the same argument could be repeated in presence of other times 
tj, giving a priori more restrictions on the position of the singularities. Eventually, the 
singularities of Wn can be inferred from their definition (Eqn. 1-17) and the singularities 
of ip'^^\ We just notice that this definition does not introduce poles at coinciding points. 

2.4 Conclusion 

We may formulate a stronger version of Theorem 1.1. 

Theorem 2.1 Assume that: 

(i) L depends on some parameter N, and has a limit when N oo. 

(a) The spectral curve of the system Eqn. 1-12 has a large N limit Sqo which is 
regular, and has genus 0. 

[Hi]' Wi(x) admits an asymptotic expansion when N oo of the form : 

Wi = Tligxi^^""^^^!! ^''^d Wf{x) for g may have singularities only at 
branchpoints of Sqo ■ 
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Then, VV„(xi, . . . , Xn) admits an expansion of the form: 

Wn = Y, iV'"''"" (2-26) 

The expansion coefficients of the correlators have only singularities at the branchpoints 
of Soo (for 2g — 2 + n > 0), and are computed by the topological recursion applied to 

W^{x{z,), x{zn))dx{z,) ■ ■ ■ dxizn) = coi{J:^){zi, . . . , z^) (2-27) 
In practice, hypothesis {Hi)' requires to check that: 

• Wi has an expansion in odd powers of A^. We think that it can always be achieved 
in an appropriate rescahng of the parameters of L with A^. 

• Wf has only singularities at branchpoints of Soo- This is not too difficult, because 
we know where singularities can arise in \l/ (Section 2.3.2). 

A similar result holds when Sqo is not of genus under an extra hypothesis on Wf , 
but is beyond of the scope of the present article. 

3 Application to the Tracy- Widom GUE law 
3.1 Tracy- Widom GUE and integrability 

The Tracy- Widom GUE law Fgue(s) can be defined as a Fredholm determinant: 

FGUE(5) = det(l-KAi)^.([,,^[) (3-1) 
where Kai is the Airy kernel: 

Jk x-y 

In their celebrated article [32], Tracy and Widom have proved an alternative formula 
making the link with an integrable system: 



FGUE(s)=exp(^y H{t)dtj (3-3) 

where H{s) = —q'^{s), and q{s) is the unique solution of the Painleve II equation: 

q"{s) = 2g3(s) + sq{s) (3-4) 

which satisfies [25]: 



^(^) - Ai(.) r. ^fy^ (3-5) 
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H{s) can be identified with a Hamiltonian for PII [29], and Fgue(s) with a r-function 
associated to this family of Hamiltonians ([24] in the sense of Okamoto, [6] in the sense 
of Jimbo-Miwa-Ueno). 

The Painleve II equation appears [22] as the compatibility condition of the following 
system for \E'(x, s): 

d^^ = 

ds^ = (3-6) 
= \l/ exp [i (|a;^ + sx) as] when x — > +oo 

where \E' = 1 + 0{l/x) when x — )■ +oo. The Lax pair is given (L, M) is given by: 

, X _ ( -4ix^ - i{s + 2g2(s)) 4xg(s) + 2ip{s) 

~ V ^xq{s) -2ip{s) 4ix2 + i(s + 2g2(s)) 



q{s) ix 

The necessary condition of existence of is S^L — S^jM + [L, M] = 0. This implies 
that q{s) is solution of PII, p{s) = q'{s). The asymptotic behavior of \& determines 
the asymptotic behavior of Eqn. 3-6 for q{s), picking up the Hastings-McLeod solution 
of PII [25]. The existence of the solution \E'(x,s) was shown in [22]. This system 
is an example of isomonodromy problem, itself closely related to a Riemann-Hilbert 
problem. Many authors have contributed to the study of these systems, and we refer 
to [23] for a review of the theory. 

3.2 s — > — oo asymptotics and spectral curve 

Let us introduce a redundant parameter A^. We define: 

q{s) = N'/'Q{S) (3-8) 



X = N^I^X 

S = iV2/3 S 



Then, Eqn. 3-6 is equivalent to: 




with the Lax pair: 
L(X,^) 



(3-9) 

$ exp [iN (|X3 + SX) as] when X ^ +oo 



2iQ'{S) 



-4tX^ - i{S + 2Q\S)) 4XQ{S) + 
AXQ{S) - AtX^ + i{S + 2Q\S)) 



MiX,S) ^ (-^^ ) (3-10) 
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and the compatibility equation: 



1 

iV2 



Q"{S) = 2Q{S)' + SQ{S) 



(3-11) 



We are now in the situation described in Section 1, where each derivative appears with 
a prefactor 1/A^. Again, Q{S) is given by the Hastings-McLeod solution to Painleve II, 
which is the unique solution [25] such that, for 5* < 0, limAr_i.oo Q{S)'^ = — f • Because 
of this property, we may study the s — ?■ — oo asymptotics for the system Eqn. 3-6 by 
studying the N —> oo asymptotics of the system of Eqn. 3-9. Then, the basic remark is 
that the derivative terms are subleading. After Eqn. 1-20, the finite spectral curve 
for this system is: 



y2 = -16X^ - SX'^S -S^- AQ\S) - ASQ^S) + 



(3-12) 



ATS 



In the large A^ limit, since Q^{S) 



S/2 (we assume S < 0), we obtain: 




(3-13) 



It can be brought in a canonical form by rescaling: 



(S) : f = hi;^{x' + 4) 



(3-14) 



where 




For this transformation, YdX = (— S'/2)^/^ ydx. Going from the topological recursion 
of Soo to that of S is only a matter of rescaling: 



F5(Soo) = (-S/2)3(1"^')F5(S) (3-16) 
4^)(Soo) = (-5/2)3(i-^-"/2)a;;i(S) (3-17) 



Thus: 



= ^(Ar2/3(_5/2))3(i-^)F3(S) 




(3-18) 
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3.3 Existence of the 1/N expansion 

In our case, L(X, 5*) has only one pole, at X = oo, which is already present in 
L*^'^)(X, S). We define Y{x) by Eqn. 3-13 with the choice of the branch of the square 
root imposed by the X — )■ oo asymptotics in Eqn. 3-9: 



r(X) = -4iX WX2 + ^, Y{x) ~ -4iX2 (3-19) 
V 2 X— ^+oo 

In order to identify J-'(Soo) computed in Eqn. 3-18 with the s — )■ —oo asymptotic of 
the T function of the initial system (Eqn 3-6), we have to check that the assumptions 
announced in Sections 1.3-1.5 hold. 

Remark 3.1 Sqo is a regular genus spectral curve: it has two simple branchpoints, 



at X = ±a/— 5/2. It admits a rational parametrization: 

X{z) = 7 + J) , y{^) = f - ^ ) (3-20) 

where 7 = \J~^- Notice that it is not the spectral curve of a matrix model, since there 
is no 1/ z term in Y . This is rather the curve of the limit of a matrix model. 

Lemma 3.1 '?/'(X, S*) (resp. ip{X,S), (f>{X,S), (p{X,S)) admits a 1/N expansion: 

'4 



^(X,5)= ^^(°)(X,5) + J^X-'V^(')(X,5)j exp 



iN ( ^X^ + SX 



(3-21) 



and for all I > 1, ^/'*^'^(X, 5*) has only poles at X = ±y (the branchpoints of Sooj. 
In particular, it does not have singularities at X = 0, the other zero ofY. 

proof: 

We already know that Q^^\S) = \J~^i and Eqn. 3-11 implies that Q has a l/X^ 
expansion: 

Q{S) = g(°) (5) + N~^' Q« {S) (3-22) 

Z>1 

Hence, L and M have a 1/X expansion. Let us apply the discussion of Section 2.3, 
where S plays the role of the time t. The large X limit spectral curve Sqo associated 
to the system dx"^ = L\E' is: 



Y{X) = ±4zXVX2 + S/2 (3-23) 
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and ^p to leading order is given by: 



cte(°)(^) 



-5/2 



X2 + S/2 



exp 

1/4 



X 



exp 



(3-24) 



To agree with the asymptotic of Eqn. 3-9, we must choose the minus sign. On the other 
hand, the large N hmit spectral curve Tf_^ associated to the system ds"^ = M\E' is: 

Y{X) = ±V-detM(o) = ±i^X^ + S/2 (3-25) 

Hence, Y{X) and Y^{X) have no common zeroes (apart from the branchpoints X = 
±a/— 5/2). So, to all orders in 1/A^, ip cannot have singularities at X = 0. The same 
discussion would hold for 0, and 0. □ 



Corollary 3.1 VV„(Xi, . . . ,X„) admits a l/N"^ expansion: 

>V„(Xi, . . . , X„) = 5^ N^-^^-^ >V^(Xi, . . . , X„) 

3>0 



(3-26) 



and the singularities of W^(Xi, . . . , X„) ati;ay from Xj = oo are found only at branch- 
points Xj = a/— 5/2. 

proof: 

There exists at least an expansion in 1/X, and the position of the singularities at all 
orders is a consequence of Lemma 3.1. We have seen in Thm. 2.1 that W„ starts 
as 0(X^~"). It remains to prove that VV„ has parity (— 1)*^ in X. Let us stress the 
dependence in X by writing Ljy, and n for the solution of Eqn. 3-9. We observe 
that *L_Ar = Ltv, which implies that *^E'l]v obey the same differential system as ^. 
Moreover, *\E'ljv has the same asymptotic behavior near x — )■ oo as \E'Ar, and is also of 

^^AT, and at the level of the integrable kernel: 



determinant 1. So, *\l/_]v 



Xi - X2 

0_jv(Xi)^/'_jv(x2) - (I)-n{Xi)^_n{x2) 



Xi - X2 



}C_n{x2,Xi] 



But, we can see on definition (Eqn. 1-17) that the correlators take a (—1)" sign if we 
revert the orientation of all the cycles. Thus: 



Wn{Xi, Xn)-N = (-l)"W„(Xi, . . . , Xn)N 



(3-27) 
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□ 

Accordingly, we can apply Thm. 1.1, which says: 
Wl^\X{z{),...,X{z^))dx{z^)---dx{z2) = a;^^(S„o)(^i,...,^„) 

(3-28) 

3.4 Tau-function and symplectic invariants 

The next step is the computation of the r-function. Let us do the computation directly, 
as an illustration of the general proof given in Section 1.5. The r-function associated 
to the unique solution of the system satisfying Eqn. 3-9 has: 

Too (x, s) = iN + SX^ a3 (3-29) 

This is correct for any A^, we see on this example that Too is indeed given by its large 
N limit. The r-function of Jimbo-Miwa-Ueno is given by (Eqn. 1-39): 

9,, In r = 2tN Res dX X Wi (X) (3-30) 

X— 5-00 

Let us compare with the variation of J-'(Soo) wrt S, given by Eqn. 1-3. We have to 
represent dsYdX with the Bergman kernel: 

B{z„z,) = y^^^ (3-31) 



We write: 



{dsYdX){z) = ^;£ffi££a = -zdJ Jx^(z) + ^ 
^ ^ y/X^{z) + S/2 I V ^ ^ 2 



Thus: 



-t Res B{z, WX2(C) + ^ = -Z7 Res B{z, C) (C - 7 
= ^7 Res B{z,C) f C - 7) 

dsF^{E^) = Z7 Res (c - I) l^f (X(C))dX(C) (3-32) 

Soo is an hyperelliptic curve, with involution C i-> ^. From construction in the topo- 
logical recursion, we know that: 

Wf (X(C)) + Wf (X(l/C)) = (2y) + (X(C)) (3-33) 
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where y+ is the divergent part of Y{X) when X — )■ +00: 



c + 

2 V C 

With the last expression, one may check that adding to Wf in Eqn. 3-32 do not 
change the result. In a first step, we may replace Wf by: 

Wf = Wf - 6g,oy+ (3-34) 

in Eqn. 3-32. Now, we have the symmetry VVf — )■ — VVf when ( — )■ This implies 
that the residues at and 00 are equal: 

dsF>^{j:^) = 2^7 Res (( - l) Wf (X(C))dX(C) (3-35) 



In a second step, we can go back to Wf : 

dsF%E^) = 2z7 Res (( - Wf{E^){X{0)dX{() (3-36) 

Another property of Wf from construction is that it behaves as 0{l/() (and even 

0(1/C^) for g > 1) when ^ — )• 00. Accordingly, we can replace — by (^C + ^j- 
Eventually: 

dsF^iJ:^) = 2^7Res fc+7Vf(Soo)(X(C))dX(C) 
= 2i Res dXXWf(Soo)X 

And, for the full series: 

a5^(Soo) = 2zX Res dX X Wi(X) = 5s In r (3-37) 



X— s>oo 



3.5 Relation with the 1-matrix model with hard edge 

So far, we have proved that the Tracy- Widom GUE law Fgue(s), defined axiomatically 
from the Airy kernel (Eqn. 3-1), has the following s — )■ —00 expansion: 

Fgue(s) = C exp I ^(-s/2)3(i-^) j (3_3g) 

\g>0 J 

where C is a constant, and S is the spectral curve of equation y'^ = |X^(x^ + 4). It 
admits the rational parametrization: 
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To complete the definition of tlie spectral curve, let us say that the physical sheet is 
associated to x = +00 [z — ?■ —zoo), and that we take the Bergman kernel: 

In [7], starting from a matrix model with quadratic potential, where all eigenvalues 
are restricted to be smaller than a given a, we obtained heuristically the following 
s — — 00 expansion: 

Fgue(s) = Ctw exp ( ^(-s/2)3(i-^)F3(Stw) ) (3-41) 

\9>0 / 

with a computable constant Ctw = 2^/^^e^'(~^'', and Stw is the spectral curve of 
equation y'^ = x + - — 2. It admits the rational parametrization: 



~ c 

We associate the physical sheet to C +00, and we take as Bergman kernel: 



5tw(Ci,C2) = 7^^^ (3-43) 

(Ci - (2)^ 

Actually, these two spectral curves are related by a symplectic transformations, 
under which the F^'s are invariant, as explained in Fig. 1. Thus, we have justified 
Eqn. 3-41, as announced in Prop 0.1. 



{x y) 



<0 = (:-\ 



reparametrization 
y{z) = -z^ " 



x[z) 



Figure 1: Symplectic equivalence of Stw and S 
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3.6 Remark on correlation functions 

We insist on the fact that symplectic transformations do not conserve the correlation 
functions. A priori: 

cu^KEtw) + ujlij:) (3-44) 

Nevertheless, one knows [19] that their difference is an exact form in each variable. 
It was argued in [7] that: 

c^n(STw)(Ci, • • • , Cn) = J](_s/2)3(i-3-"/2) a;^^(STw)(Ci, • • • , Cn) (3-45) 

is the asymptotic expansion when s — )■ — oo of the (rescaled) correlation function of 
eigenvalues \i = 2 + sN^'^/'^{l -l- y). In other words, for s < and any Borel subsets 
Ji, ■ ■ ■ , Jn ] — oo, 0[, we have: 



lim Prob [Ai e Ji, . . . , A„ G J„ I A^ax < 2 + sA^-^/si 

Wn(STw)(Cl,---,Cn) (3-46) 



iV-s-oo 

dCi / dC 



in the sense of an asymptotic series when s — t- — oo. Ki is the image on the curve of Jj. 

A priori, this is different from the determinantal point process defined with the 
integrable kernel: 

^ ^(^1)^(^2) -^(xO0(x,) 

of the system of Eqn. 3-6. 
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4 Conclusion 



Tracy- Widom law 

Let us summarize what we have obtained. 

• Near an hard edge a approaching the natural soft edge amax, the density of 
eigenvalues of a random hermitian matrix is (after appropriate rescaling) locally 
given by 

(Stw) : y=\ x + -- 2 
V X 

and therefefore, the partition function in this regime is asymptotically given by 
the symplectic invariants of Stw- 

InZ(a) ~ ^(1 - a/a^a.)(2-2^)i N'-'^ F^^tw) 

9 

and the eigenvalues correlation functions of the random hermitian matrix are 
also asymptotically given in this regime by the symplectic invariant correlators 
of Stw, i-e. by the uJn\'ETw)- Let us also mention that all this holds for any 
problem, solvable by symplectic invariants, whose spectral curve is locally given 
by Etw- This is for instance the case for the statistics of plane partitions near 
the boundary of the liquid region [21]. 

• The second step is that Etw is symplectically equivalent to the spectral curve S: 



-xy/ x'^ + 4, 

which shows that 

InZ(a) ~ - a/a^ax)^'"'^^^ N^'^^ 

9 



and since S is the classical spectral curve of a Lax pair for the Painleve II equa- 
tion, we recover, as expected, the large s expansion of Tracy- Widom law. 

• However, in the symplectic invariance of Sxw and S, only the F^'s are invariant, 
the correlators are not invariant. 

It is not clear to us how the symplectic transformations (in particular x i-^ y, 
which is really the non trivial one at the level of the F^) act on an integrable system of 
classical spectral curve (x, y{x)). It would be interesting to exhibit a full {N dependent) 
integrable system whose classical spectral curve is directly Stw; without having to 
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perform a symplectic transformation. Such a system should be related to PIT, because 
the T function is the same as that of Eqn. 3-6, it can be expressed in terms of a solution 
of PII. 

For (5 ensembles, the generalization of this approach is an open problem. There 
exists a version of the topological recursion and symplectic invariants for all /3 > 0, 
computing the large N expansion in the /3 random matrix ensembles [10] provided that 
the expansion exists. Yet, the analog of Section 1.3 for /3 7^ 2 is missing: starting from 
an integrable system, it is not know how to define quantities Wn''^ satisfying the loop 
equations of /3 ensembles. One could expect that (3 ensembles should be related to 
quantum integrable systems, having to do with Bethe ansatz, instead of classical Lax 
pair systems. 

Loop equations and integrable systems 

We have reviewed that it is possible to associate loop equations (or Virasoro constraints) 
to any first order 2x2 rational differential integrable system. In fact, we have completed 
the work of [3] by: 

• Showing that hypothesis on analytical properties and 1/A^ expansion of the n- 
point functions need only to be checked for the 1-point function. 

• Clarifying when these analytical properties and expansion properties holds. 

• Under all these assumptions, providing a proof of the former claim that the 
summation of symplectic invariants reconstruct the Tau function. 

To fulfill this program, it was essential to include the differential system in an isomon- 
odromic problem with a full set of times. It is in fact possible to associate loop equations 
to first order system of size d x d, and we are currently working on the generalization 
of the same theory to these systems [8]. 
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A A counterexample to the expansion 

Consider a system 1/N d^'^ = L\l/. The first loop equation is: 

W2{x,x) + (Wi{x)f = N^Pi{x) (1-1) 
where Pi (x) is a rational fraction given by: 

Pi = -detL = ^TrL^ (1-2) 
Let us assume that L = L(°) + ^L^^) + 0{1/N^). Then: 

Pi = ^TrL(°)' + i^TrL(°)L(i) + 0(1/Ar2) ^^,3) 

We want to find a case where the W„ do not have an expansion with fixed parity of 
1/A^. It is enough to exhibit a system for which the first two terms in Pi do not vanish. 
We call Painleve IIq, the equation: 

q"{s) = 2q{sf + sq{s) - a (1-4) 

where a is a fixed parameter. This equation appears as the compatibility condition of 
the Lax system [22]: 

(9^* = 



(9,* = 
Lfx, s) = 



s 

2 „Vo I 0^2 



M{x,s) 

Let us do the rescaling: 



-Aix^ - i{s + 2q^{s)) Axq{s) + 2iq'{s) + f 

Axq{s) - 2iq'{s) + f 4ix^ + i{s + 2q^{s)) 

—ix q{s) 
q{s) ix 



,Jn2/3s ' q{s) = N'/'QiS) (1-5) 



This leads us to the compatibility equation: 



1 rv 

— Q"{S) = 2Q{Sf + SQ{S)-- (1-6) 



and the Lax system: 



N 



, -4tX' - ^{S + Q\S)) AXQ{S) + m^ + j^ 
^ ' ' ' 4Xg(5)-^ + ^ AtX^ + tiS + 2Q\S)) 



M{X, S) 



-iX Q{S) 
Q{S) iX 
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When — )■ oo, we have to the first two orders 

1 a 



Q{S) = y/-S/2 + ^ ^ + 0{1/N') (1-7) 



We may compute: 
TrfL^L^i) 
= Tr 

7^ 



Hence, the fixed parity property is not satisfied. Yet, one could also make the choice to 
rescale a in A = a/N, and keep A fixed. In this case, we could reproduce the argument 
of Section 3.3, and check the fixed parity property. This counterexample supports the 
idea that a l/N"^ expansion may exist only in a good choice of rescaling with N of the 
parameters of the differential system. 

B Topological recursion and symplectic invariants 

With the data of a spectral curve (see Section 1.1), the topological recursion defines a 
tower of differential forms as follows: 

Definition B.l We define: 

oj['\z) = -y{z) dx{z) (2-1) 

J^\zi,Z2) = B{zi,Z2) (2-2) 

and by recursion on 2g — 2 + n: 

J 



i 

+ E E ^) J \ I)] (2-3) 

h=Q ICJ 

where the sum over i is over all branchpoints ai (the zeroes of dx{z)), residues are 
taken as contour integrals on the Riemann surface C, the kernel K{zq, z) is: 

KM^~ (2-4, 
2(4"'(E)W-c.<°'(E)(z)) 

and in the sum Yl'icJ prime ' means that we exclude the pairs (/i = 0, / = 0) 
and {h = g,I = J) from the sum. 
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An important property is that for 2g — 2 + n > 0, Un'^ is a meromorphic form in each 
Zj, with poles only at branchpoints, and is symmetric in all Zj's. For instance, we have: 

uj'i-\zQ) = Res K{zo,z) B{z,z) (2-5) 

i 

col'\zo, zi, Z2) = y2 Res Kizo,z){B{z,zi)B{z,Z2) + B{z,Z2)B{z,zi)) (2-6) 



Then, having defined lJ£^ for n > 1: 



Definition B.2 For g > 2, we define the symplectic invariant uJq^\t,) also denoted 
F^'(S) as: 

F^(S) = 4^HS) = -i- ^ Res Hz)Jf\E){z) (2-7) 

^ i 

where ^ 

<l(z) = / y{z')dx{z') (2-8) 



and the value of is independent of the base point a and path used to compute the 
integral. 

The definitions of F°(S) and -F^(S) are more involved and we refer the reader to [19]. 
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